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Motivated by a recent experiment by Zelevinsky et al. [Phys. Rev. Lett. 96, 203201], we present the theory 
for photoassociation and optical Feshbach resonances of atoms confined in a tight one-dimensional (ID) or 
two-dimensional (2D) optical lattice. In the case of an alkaline-earth intercombination resonance, the narrow 
natural width of the line makes it possible to observe clear manifestations of the dimensionality, as well as 
some sensitivity to the scattering length of the atoms. Among possible applications, a 2D lattice may be used to 
increase the spectroscopic resolution by about one order of magnitude. Furthermore, a ID lattice induces a shift 
which provides a new way of determining the strength of a resonance by spectroscopic measurements. 



I. INTRODUCTION 

There is a growing experimental effort to study the proper- 
ties of ultracold alkaline-earth vapours. One of the reasons is 
that the narrow intercombination resonance of alkaline-earth 
species (weakly coupling 1 S'o and 3 Pi states) can be used 
to build optical atomic clocks which could be more accu- 
rate than the current atomic standard of time JH 01 • The 
Bose-condensation of Ytterbium |4], whose atomic structure 
is close to alkaline-earth species, has also raised the hope 
of condensing these species. To reach these goals, a good 
knowledge, and possibly control, of alkaline-earth atomic in- 
teractions are needed. Photoassociation, the process of asso- 
ciating pairs of colliding atoms into excited bound states by 
making them absorb a resonant photon, appears as the best 
tool to characterize and control these interactions. Indeed, 
photoassociation can be used as a spectroscopic tool to mea- 
sure energy levels in excited molecular states 1 5J |fj [J and 
characterize ground-state atomic interactions |8]. It can also 
be regarded as an optical Feshbach resonance |9], analogous 
to magnetic Feshbach resonances in alkali systems, making 
it possible to alter these interactions I10H . This is particularly 
useful for alkaline-earth systems with no hyperfine structure, 
where magnetic Feshbach resonance is not possible. 

Photoassociation near the intercombination resonance is 
characterized by narrow lines which are sensitive to effects 
like recoil shifts (due to the momentum of the absorbed pho- 
ton) and Doppler broadening. These effects can be elim- 
inated by strongly confining the atoms in the direction of 
propagation of the photoassociation light. This has been 
demonstrated by trapping the atoms in tight optical lattices 
The wave length of the lattice can be ad- 
justed to its "magic" value lUTll . so that atoms in the ground 
state and the excited state feel the same trapping potential. 

In turn, it is known that strong confinement may affect the 
collisional properties of the atoms IU41 11511 . We can there- 
fore expect photoassociation to be different in an optical lat- 
tice than it is in free space. The purpose of this paper is 
to investigate how photoassociation spectroscopy and opti- 
cal Feshbach resonances are affected by a one-dimensional 
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(ID) or two-dimensional (2D) optical lattice. In Section ITU 
we present the general theory of resonant collisions in a ID 
and 2D optical lattice, building upon the works of D. S. 
Petrov and G. V. Shlyapnikov Hj and M. Olshanii 0- 
In section we apply the theory to the case of alkaline- 
earth species. In section ffvl we conclude on the possible 
effects and uses of alkaline-earth photoassociation in optical 
lattices. 



II. OPTICAL RESONANCES IN A ID OR 2D LATTICE 

An optical lattice is obtained using a standing wave laser 
which creates a sinusoidal periodic trap. Each cell of this 
lattice confines the atoms in either a pancake-shaped cloud 
(for ID lattices) or cigar-shaped cloud (for 2D lattices). For 
this reason, we will refer to ID lattices as a 2D (pancake) 
geometry, and to 2D lattices as ID (cigar) geometry. We will 
assume that there is little transfer between each cloud, so 
that we can treat each one independently. The atoms in each 
cloud will be assumed free to move in the directions which 
are not confined by the lattice, ie, the transverse directions 
x, y for a ID lattice, and the axial direction z for a 2D lattice. 
We will also describe the confinement induced by the lattice 
near the centre of each cell by a harmonic potential, which is 
valid if the lattice is tight or strong enough. 

To describe optical resonances under such confinement, 
we use the theories of collisions in ID or 2D geometries of 
Refs. We express these theories in paragraph !!! Al in 

terms of an energy-dependent scattering length a(k), which 
is essentially the i^-matrix for 3D collisions. The origi- 
nal theories were formulated in terms of the usual (energy- 
independent) scattering length, but taking into account the 
energy-dependence extends the range of validity of the theo- 
ries I16T . and is necessary in the case of resonant collisions. 
We then give the expression for a(k) in the case of resonant 
collision in paragraph 111 Bl Using this expression, we de- 
duce the collisional and photoassociation rates in ID and 2D 
geometries. 

Similar methods have been used to treat specific cases: a 
A"-matrix approach was used in Ref. 1 1 711 to describe bosonic 
resonant collisions in a 2D geometry, a multi-channel K- 
matrix approach was used in Ref. 1 18] to describe collisions 
of spin-polarized fermions in a ID geometry, and a renor- 
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malization approach was used in Refs. to describe 

resonant collisions of bosons in a ID geometry. All these 
works relate the effective ID or 2D scattering properties to 
the free-space properties. Also closely related are the con- 
structions of pseudo-potentials which can treat confined sys- 
tems Il2ll,l22l 12311 and pseudo-potentials in low-dimensional 
systems I22M24TI . 



energy in the weak direction: 
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A. Collisions in a ID or 2D confined geometry 

Let us briefly review here the method of Refs. fl4l fl^l . 
The stationary states i/i(r) for the relative motion of two 
atoms confined by a harmonic potential V(r) are given by 
the Schrodinger equation 



-—V 2 + U(r) + V(r) 



ip(r) — Eip(r) 



(1) 



where r = (x,y,z) is the relative coordinate between the 
two atoms with length r = \r\, fj, is their reduced mass, U(r) 
is their interaction potential, E is their relative motion en- 
ergy. In a 2D geometry (ID confinement) the trapping po- 
tential is: 

V{r) = V 2D (r) = ^n 2 z 2 

where f2 is the frequency of the trap. In a ID geometry (2D 
confinement), it is: 

V{r) = V lD {r) = \^ 2 p 2 

withp = (x,y). 

At separation r larger than the range r$ of the interaction 
potential U, the stationary states can be written as a linear 
combination of products of a harmonic oscillator eigenstate 
in the confined direction and a free wave in the weak direc- 
tion. We choose the conventional basis where each station- 
ary state in the weak direction is asymptotically the sum of 
an incident plane wave and a scattered wave (which is as- 
sumed isotropic for low energies). In 2D and ID, we have 
respectively: 
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where ip v , cf> nm are respectively the unit-normalized ID (2D) 
harmonic oscillator wave function of vibration index v (prin- 
cipal number n and angular momentum number to), p is the 
wave number of the incoming ID wave, and if is the wave 
vector of the incoming 2D plane wave with wave number q. 
For each state, the total energy E is composed of the oscil- 
lator energy in the confined direction, and the free particle 
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which defines the wave numbers p and q for the incoming 
plane wave, as well as p' and q' for the scattered wave. The 
basis of eigenstates (0 and (|2} define respectively the effec- 
tive 2D (ID) scattering amplitude f^ v , Unm,n'm'^ for two 
atoms colliding with relative energy E in the initial state ip v 
(4>nm) to end up in the final state ip v i ((f> n 'm')- 

At short separations, one can neglect the confining trap in 
Eq. Q and the collision recovers the aspects of a usual 3D 

t-2,2 

(free space) collision at energy E = ^^j^- As a result, as 
long as kro < 1, there is a region 1 16] where the the wave 
function is proportional to the usual 3D s-wave scattering 
wave function 
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where a(k) is the e nerg y-dependent s-wave scattering length 
introduced in Ref. I25ll2al . At very low energy, the scatter- 
ing length a(0) is simply the usual scattering length a. More 
generally, it is related to the 3D elastic scattering matrix ele- 
ment S(k) by: 



»(*) = 



1 1 - S(k) 
ik 1 + S(k) 



(7) 



By matching Q or (|2} with 0, one finds the relation be- 
tween the effective ID or 2D scattering amplitude and the 
scattering length fl5H27ll : 
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where ao = y/h/fiCl is the harmonic oscillator length, w(x) 
is a function introduced in Ref. fl5ll . and £(a, x) is the Hur- 
witz zeta function (as defined in Ref. (3). These functions 
have the following properties: 



for x 1 

for 0.1 < x < 1 

for x < 1 
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for x <C 1 



where B « 0.915 and ((i) 
properties of ip v and cj> nm [l 
both and v' are even, and / 
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B. Optical Feshbach resonances in free space 

Usual photoassociation in free space can be modelled 
by a multichannel scattering theory l29Tl . In the case of a 
photoassociation resonance with an isolated molecular state, 
one finds simple expressions for the scattering properties of 
atoms under the influence of the laser light. For instance, the 
s-wave elastic scattering matrix element reads: 



S(k) = S bg (k) 1 - 



iT(k) 



E-E + UT(k)+ 7 ) 



(14) 



Here, 5b g (fc) is the s-wave elastic scattering matrix ele- 
ment in the absence of light, E = is the collision en- 
ergy, Eq is the position of the resonance (with respect to the 
ground-state threshold), T(fc) is the stimulated width of the 
resonance, and 7 is the natural width of the resonant molec- 
ular state. The stimulated width T(fc) is due to the coupling 
by the laser between the incoming scattering state and the 
resonant molecular state, and is proportional to the intensity 
of the laser light. The natural width 7 is due to the decay of 
the resonant molecular state by spontaneous emission. Be- 
cause of this natural width, the scattering matrix element is 
no longer unitary (^(/c)! 2 < 1), meaning that the collision 
is no longer elastic. This accounts for the fact that some 
photoassociated pairs of atoms can deexcite by spontaneous 
emission, and do not return to their initial state. In most 
cases, these pairs of atoms gain enough energy to leave the 
trap where the system is confined, or form bound molecules 
13011 - so that photoassociation is observed as a loss from the 
initial atomic cloud. The observed loss serves as a basis for 
photoassociative spectroscopy. 

A magnetic Feshbach resonance is described by a formally 
similar theory, where spontaneous emission is generally neg- 
ligible (7 = 0).The following theory is therefore applicable 
to a magnetic Feshbach resonance by setting 7^0 and 
7^opt — * where k is the strength of the magnetic resonance. 



From the expression dl4> of the scattering matrix element, 
one can derive the energy-dependent scattering length Q 
near a photoassociation resonance: 
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where a Dg (k) — i(l — Sb g (fc))/(1 + S ng (k)) is the back- 
ground energy-dependent scattering length in the absence of 
light, and the optical length l op t(k) is defined by [10]: 
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Wigner's threshold law 

For sufficiently small collisional energies, the incoming 
scattering state obeys Wigner's threshold law. Namely, the 
energy-dependent background scattering length ab g (fc) be- 
comes constant 



(17) 



abg(fc) ~ a bg , 

when ifc 2 ab g r e <C 1, where r e is the effective range 1 3 1 ] 
of the background interaction (usually on the order of ro). 
Here, a\, g denotes the usual zero-energy background scatter- 
ing length <2bg(0). Similarly, for kr e <C krc -C 1, where rc 
is the Condon point of the photoassociation transition, the 
stimulated width is proportional to the collisional momen- 
tum, and the optical length becomes a constant: 



^opt(^) ~ l> 



opt ■ 
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This constant l opt characterizes the strength of the resonance 
at low energy. 

In many cases, ka bg and kl opt are small, so that Eq. d!5l > 
simplifies to: 



a(k) « a bg + Z op t 
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which is the usual expression for the scattering length of the 
atoms modified by the radiative coupling. 



C. Resonance in a 3D geometry 

In a 3D (almost free space) geometry, the number of (col- 
lision or photoassociation) events per unit of time is given 
by: 



K. 3D (n 3D (x,y,z)) dxdydz 



(20) 



where n 3D is the density of atoms and K, 3D is the usual 3D 
(collision or photoassociation) rate, which can be calculated 
from the scattering properties of the atoms, such as the scat- 
tering matrix element (114-1 . They are given respectively by: 
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where 



E. Resonance in a ID geometry 



.)V(k) 



denotes statistical average over all the possible collision 
wave vectors k, and'P(k) is the wave vector distribution nor- 
malised to 1 . One can see in (12 li that the photoassociation 
rate follows from the non-unitarity of the scattering matrix, 
IS^fc)! 2 < 1. From Eq. H41 . we also note that the photoas- 
sociation rate does not depend directly on the background 
scattering length <2b g . There is, of course, an indirect depen- 
dence through r or Z opt which contain a Franck-Condon fac- 
tor involving the ground-state wave function, whose nodal 
structure depends on the background scattering length. This 
indirect dependence has been used to infer the background 
scattering length from photoassociation spectra 1 8 - l32tl . 



D. Resonance in a 2D geometry 

In a 2D geometry, the number of (collision or photoasso- 
ciation) events per unit of time is given by: 



K 2D {n 2D (x,y)Ydxdy 



(22) 



where n 2D is the 2D density of atoms, and the rates of col- 
lision or photoassociation JC 2D are related to the quasi-2D 
scattering matrix element S^ v , = 5 uu i + hfy U i' 



In a ID geometry, the number of (collision or photoasso- 
ciation) events per unit of time is given by: 



ID (\D 



{z)) 2 dz 



(26) 



where n 1D is the ID density of atoms, and the rates of 
collision or photoassociation K. 1D are related to the quasi- 
1D scattering matrix element S® mn , m , = 8 nn >5 mm , + 

J n'm,n'm' 



'Vol 



2/i 




\ n m i n,m,p 

(27) 

In these expressions, a sum is performed over all final 
states nf, to', and the average 



n,m,p) 



denotes the statistical average over all initial states of relative 
motion. The statistical distribution V(n, m,p) is normalised 

to unity: fdp Y, n ,m P ( n > m 'P) = L 

From Eq. (|9jl and the properties of £, one can derive: 
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In these expressions, a sum is performed over all final 
states v' , and the average 



d 2 q 



denotes the statistical average over all initial states of relative 
motion. The statistical distribution V(v, q) is normalised to 
unity: Jd 2 q "Z„V(v,q) = 1. 

From Eq. and the properties of w, one can derive: 
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where we recall that k 2 = 2j ^- 

a 

factorial of n. 



q 2 and n!! is the double 
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F. Comparison of 2D/1D and 3D theories at thermal 
equilibrium 

The comparison of 2D or ID rates with the more familiar 
3D rates is not direct. It is more easily achieved at thermal 
equilibrium. First, we can evaluate these rates assuming a 
Boltzmann distribution in all cases: 



V(k) 



a 3 / 2 exp 
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where a = 2 -K^,k B T anc ^ P = 1 — e k s T . Secondly, in a 
strongly confined gas at temperature T, the density profile in 
the confined direction can be assumed to be "frozen" in the 
Boltzman profile (at least, for a certain time): 

n 3D {x,y,z,t) « n 2D (x,y,t) x f T (z) 
n 3D (x,y,z,t) « n 1D (z,t)x f T (x)f T (y) 

with the Boltzmann profile: 

f T (x) = —=-e-^y 



where L — a / ^/2 tanh(M7/2fc B T). Note that a thermal 
distribution of independent distinguishable atoms is assumed 
here. Then we can state that the number of events per unit 
of time (12 2 1 or i26\ is given by the 3D formula J20i with an 
effective 3D rate JC 2D ', or K, 1D respectively. Integrating out 
the profile in the confined direction, we can find the relation 
between the 2D/1D rate and the effective 3D rate: 

K? D = V^L xJC 2D , 

IC 1D = (V2Hl)\jc 1d . 

The effective 3D rate thus calculated from the 2D or ID the- 
ory can be compared with the usual 3D rate Oil . In the 
following, we consider three regimes identified in Ref. fiUl . 
corresponding to the different regions \\2\ of the function w 
or £. 



Free space regime ksT ^> hfl 



Confinement-dominated regime O.lftfi < ksT < hQ. 

In this regime, very few incoming states in the confined 
direction contribute to the collision. We can make the ap- 
proximation that only the v = 0, or n, m = incoming 
states make a significant contribution. In a 2D geometry, the 
effective 3D rates become: 




where ( • • • ) — / d 2 q (. . . ) V(q] is the statistical aver- 
age over all 2D momenta, V(q) = 2nuk^T ex P ( ~ ksT ) 
is the 2D Boltzmann distribution, k = J a^ 2 + q 2 , and 

A = y 1 — e fc B T isa correction factor close to 1. Com- 
parison with ( 13013 1> shows that the ID confinement brings 
two main effects: 

• the 3D momentum integral becomes a 2D integral 

• the momenta involved in the integral now have a lower 
bound given by the "zero-point momentum" K/ao. 

Note that k is replaced by J a^ 2 + q 2 in a(k), and 
by y/ir/ao elsewhere. 



In this limit, a large number of incoming v, or n, m are in- 
volved: the states of relative motion in the confined direction 
form a quasi-continuum as in free space. Thus, we expect to 
retrieve the free-space results. Indeed, both in the 2D and ID 
theory, the effective 3D rates simplify to: 
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(30) 
(31) 



One can check that these expressions are exactly the same as 
the usual 3D rates ( 12 lb . In particular, if the Wigner's thresh- 
old laws (1171 1 8i apply (ie if the temperature is not too high), 
the photoassociation rate is given by the standard expression: 



JC, 
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where E — The term kl opt is responsible for the so- 

called power-broadening by the laser. 



In a ID geometry, the effective 3D rates become: 
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where 



x Im [— a(k)] 
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(36) 



v 



= I-oodp(- ■ -)V{p) is the sta- 



tistical average over all ID momenta, V(p) 



1/2 

2-Kfik B T I eX P 



k B T 



is the ID Boltzmann dis- 



2a n 2 + p 2 . Comparison with ( 13013 H 



tribution, and k — v 
shows that the 2D confinement brings the following effects: 

• the 3D momentum integral becomes a ID integral 

• the momenta involved in the integral are now typi- 
cally larger than the "zero-point momentum" ^/2h/ao- 

■ p 2 in a(k), and by 



Note that k is replaced by y 2a 2 
2 /(pa 2 ,) elsewhere. 
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• there is an extra term which is real and proportional to 
the Hurwitz zeta function. 

Both in ID and 2D geometries, it is interesting to note that 
the photoassociation rate now depends directly on the back- 
ground scattering length a\, g through a(k) - see Eq. ( 1151 . 



unlike in the free space limit. This dependence involves the 
ratio <2b g / ao of the background scattering length over the har- 
monic length. In many cases, this ratio is small and can be 
neglected, so that the Wigner limit ( 11711 8t of Eqs. J34l > and 
d36l simplifies to: 



hi 



x e 



dE 
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The result d37i was stated in Ref. fl3ll . omitting the factor 
A and the power-broadening term y^opt/ao- It turns out that 
Eq. ( I37t can be integrated analytically: 



£2D =A ^opt / 7 



Fn 
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where F(x,y) — ^Im [e ^ i!/ )(Ei(ai — iy) + in)] and 
Ei(z) = — J^ie^ 1 /t)dt is the exponential integral defined 
with a branch cut discontinuity from — oo to 0. 



Quasi-2D and quasi-lD regime Ml 3> fcsT 

At very low temperature fc^T < Ml in the 2D geometry, 
we need to include the logarithmic term coming from w - see 
Eqs . (I12> - which gives the character of a true 2D collision 
11511. This leads to: 
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As explained in Refs. fl5l l33ll . the logarithmic term leads 
to a "confinement-induced" resonance when it cancels the 
term ao/a(l/ao), which can be interpreted as the appearance 
of a 2D bound state at threshold. Since we assumed that 
a(k) already has a resonant character - Eq. ( H5L we may 
reinterpret the logarithmic term as a shift in the position of 
this resonance. Note that Eq. (1391 was already considered 
in Ref. fPUl to describe a magnetic Feshbach resonance in a 
quasi-two-dimensional gas. 



At very low temperature in the ID geometry, the zeta term 
in Eqs. (I35l36i becomes close to C(i) : 
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As explained in Refs. 

GUSH, the term £(§) leads to 
a "confinement-induced" resonance when it cancels the term 
ao / a(y/2/ ao) . This again can be interpreted as a shift in the 
position of the optical resonance. 

While the logarithmic and zeta terms in Eqs. J40i and ( I42t 
cause a shift in the position of the resonance, the imaginary 
terms i x /n and i- 2 - cause a broadening of the resonance 

V pan & 

which reduces the maximal photoassociation rate. Note that 
in a quasi-2D geometry, the resonance is increasingly shifted 
by the logarithmic term when the temperature is decreased, 
and it is broadened by a constant term iy/n. In contrast, in 
a quasi- ID geometry, the resonance is shifted by a constant 
term C( | )> but is increasingly broadened by the term i-^ as 
the temperature is decreased. 

The reduction of the photoassociation rate caused by this 
strong broadening in ID is related to the "fermionization" 
of bosons fl4l 13511 at low temperature, and was recently ob- 
served experimentally l36ll . This can be understood as fol- 
lows. At very low temperature the ID photoassociation rate 
Eq. E3 can be written as 



JCiP = 



-pa 



= (l^| 2 ) 



3D 



X ^na 



(43) 



where K.^P 
pa 



2h 



Im(- 



o) is the 3D photoassociation rate 
Eq. J3 11 in the limit T — > 0, and r\\ B is the renormalisa- 
tion factor ( II \\ of the relative wave function at short sep- 
aration r. The average (|7yi£)| 2 ) therefore corresponds to 
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the pair correlation function g 2 (of the two-body problem) at 
short separation. Fermionization occurs when this quantity 
becomes small, ie the probability of finding atoms at short 
separation is strongly decreased. As photoassociation typi- 
cally occurs at these short separations, the photoassociation 
rate J43l > is strongly reduced. For the Boltzmann distribu- 
tion V(p) that we have assumed so far, {\rj\D\ 2 ) p actually 
vanishes at T = 0. A more proper calculation using the 
relative momentum distribution 113 711 in the fermionization 

regime leads to {\iji D \ 2 ) p = |tt 2 (j^s) > where n iD is 
the ID density of the gas, in accordance with the calculation 
of the pair correlation function g 2 in the many-body frame- 
work of Ref. 13811 (note however that a denotes the modified 
scattering length ( 119b . instead of the background scattering 
length). 



G. Position of the resonance 

In a typical experiment, the collision or photoassociation 
rates are measured as a function of the laser frequency. In 
a lattice, the frequency for which the laser is resonant (ie 
which gives a peak in the rate) may be shifted from the ex- 
pected free space resonant frequency. The lattice eliminates 
the recoil shift but also introduces new shifts. In our equa- 
tions, the free-space resonance is obtained for Eq — 0, there- 
fore we define the lattice shift AE as the value of Eq which 
maximizes d25l > or \29\ . 

There are several contributions to the shift. The first one 
is the zero-point energy, which can be easily understood as 
follows. The relative motion of the ground-state atoms is 
influenced by the confinement, and therefore has the zero- 
point energy, even for T = 0. On the other hand, the relative 
motion in the excited molecular bound state is not affected 
by the confinement (only the centre-of-mass is). As a result, 
the transition is shifted by the zero-point energy. 

The thermal average in fl25i and fl29b also shifts the line 
by a quantity of order k B T. Due to the dimensionality of 
the integral involved in the averaging, this thermal shift will 
differ from the usual 3D thermal shift. 

In the 2D confinement-dominated regime, we find the fol- 
lowing shift: 



AE 
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CD 



-Ml + k B T x Z 



k B T/j 



(44) 



The first term is the zero-point energy shift and the second 
term is the thermal shift, where Z(y) w 1 — (2/y) 2 is the 
solution of F(Z, y) = Z i+ y 2 , and F is the function defined 
previously. 

We mentioned that the background scattering length a\, g 
affects the photoassociation rate in a lattice, and in particular 
it may shift the line if it is large enough (note that this is an 
intrinsic molecular shift, not a density-dependent shift due to 
interaction with other atoms |2]), however we find this shift 
to be very small in most cases. 

Finally, at high intensity and low temperature, the line may 
be shifted by the logarithmic term in a 2D geometry, or by 



the zeta term in a ID geometry. The 2D logarithmic shift 
is potentially the most easily observable experimentally, be- 
cause the line is not as broadened as in the ID case. We find 
that the following formula gives a good approximation of the 
shift in the quasi-2D regime at high intensity: 
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(45) 

Interestingly, the "logarithmic shift" is proportional to the 
optical length i opt . In free space, the optical length is usually 
determined by measuring the rate (1321 which is proportional 
to it. However the rate can typically have a factor of 2 un- 
certainty as its determination relies on the knowledge of the 
density of atoms M il 311 . Here we propose a determination of 
/ opt which does not require the knowledge of the density and 
is is simply based on a spectroscopic measurement: by mea- 
suring the position of the line for different intensities, one 
should observe a shift linear with the intensity I, and deduce 
the strength of the resonance l opt /I- The major constraint of 
this method is to be able to reach and measure sufficiently 
low temperatures -C M7. One should also take into account 
any possible light shift 1 39| (also proportional to the inten- 
sity) by measuring it at higher temperatures. 



H. Line shape and strength 

The shape and apparent strength of the line are also mod- 
ified by the lattice. In the confinement-dominated regime, 
we saw that the dimensionality of the thermal distribution 
is reduced. As a result, the corresponding density of states 
(~ 1 in 2D, and ~ E^ 1 !" 1 in ID) gives more weight to 
small energies - compare Eqs. d37t and (I38i with Eq. (1321 . 
The line is therefore expected to be narrower, especially in 
the ID case. The precise strength of the line depends on 
the density and the thermal averaging, but we can get some 
qualitative insight by considering the on-resonance rate in 
the confinement-dominated regime, disregarding the thermal 
average. 

In most cases, species have a small scattering length 
ctbg < In this limit, the on-resonance rates are essen- 
tially: 
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(46) 
(47) 



One can see that the power broadening term kl op t in 
Eq. d32l > is enhanced because the momentum hk is now 
bounded from below by the zero-point momentum ~ h/ao. 
It is even much larger at small temperature in the quasi- ID 
regime, because of the fermionization effect discussed ear- 
lier. It is therefore easier in a lattice to saturate the line with 
lower laser intensities, even at very small temperature. 

On the other hand, certain species may have a large scat- 
tering length dbg ~ clq. In this case, the rate explicitly de- 
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pends on the scattering length. For small laser intenstites 
Gopt <C eibg, do), the on-resonance rates are essentially: 
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(49) 



We give an illustration of this dependence on ab g in the 
next Section in the case of Strontium 86. 

In principle, this new dependence on the background scat- 
tering length could be used as a way to improve its deter- 
mination when it is large. This would bring complementary 
information to the usual node method 1 8] which proves to be 
difficult for large scattering lengths as it requires a very pre- 
cise knowledge of the long-range part of the potential. How- 
ever, the typical uncertainties associated to measured rates in 
current experiments may limit the interest of this idea. 



III. ALKALINE-EARTH PHOTOASSOCIATION IN A ID 
OR 2D LATTICE 

We illustrate the previous theory by considering alkaline- 
earth intercombinationphotoassociation in a tight optical lat- 
tice. An important feature of the intercombination line is the 
low spontaneous emission from the photoassociated state. 
While the natural width for allowed transitions is typically 
tens of MHz, the width of the intercombination lines is typ- 
ically tens of kHz. This is on the same order of magnitude 
as the typical trapping frequency at the centre of each lat- 
tice cell. As a result, we expect a clear manifestation of the 
effects described in the previous section for alkaline-earth 
intercombination transitions. On the other hand, they should 
hardly be observable with allowed transitions, as they would 
be hidden by the large natural width. 

Figure [2 shows the effective photoassociation rate at dif- 
ferent temperatures for a typical intercombination line of 
Strontium 88 in free space, ID lattice, and 2D lattice. The 
scattering length of Strontium is taken to be 10 Bohr |8] (1 
Bohr ss 0.0529177 nm) and the spontaneous emission term 
7 = 15 kHz lll3ll . The trapping frequency Q is 27rx50 kHz, 
and the optical length is set to 1 Bohr. The temperatures 
correspond to the three regimes discussed in Section 111 Fl 
free space (T = 100 /iK and T = 20 fiK), confinement- 
dominated (T = 4 [iK ), and 2D regime (T = 0.1 ). 

For T = 100 /zK and T = 20 (Fig. [TJi and[TJ»), one 
might expect the line shape to be very close to the free-space 
theory prediction. This is indeed the case in a 2D geome- 
try. However, higher temperatures are needed to reach the 
free-space regime in a ID geometry. For the temperatures 
considered here, we observe a lot of distinct spikes in the ID 
photoassociation line. These spikes correspond to contribu- 
tions from different trap states and are therefore separated by 
twice the frequency ft of the trap - see Eq. (|3J with rn = 0. 
There are similar contributions in the 2D case, but they tend 
to overlap, due to their width. In contrast, the spikes are very 
narrow in the ID case, because the one-dimensional density 



of states scales like 1/yE (see, for instance, Ea. P8l >). which 
enhances low-energy collisions. As a result each spike has a 
width close to the natural width of the resonance and is well 
separated from the others, even at relatively high tempera- 
tures (tens of fiK) which are more than one order of mag- 
nitude larger than the natural width. This feature could be 
exploited to measure the position of the line at nearly the nat- 
ural width precision, while operating at temperatures much 
larger than this width. Note that we have assumed ideal har- 
monic confinement, which is why the spikes are regularly 
spaced and well separated. Experimentally, only the first 
spikes should be visible, due to the anharmonicity of the lat- 
tice potential. 

In the confinement-dominated regime (Fig. HJ), the line 
shape is dominated by the first spike (contribution from the 
trap ground state) and the line is clearly shifted due to the 
zero-point energy. In the 2D case, the density of states leads 
to a line which is more symmetric and narrower than in 3D, 
resulting in an increase of the photoassociation rate on reso- 
nance with respect to the free-space rate. This situation cor- 
responds roughly to the experimental conditions of Ref . 11311 . 
However, the experimental conditions were such that the 2D 
and 3D theory happened to predict very similar rates, and 
extra broadening mechanisms made it difficult to observe an 
unambiguous 2D line shape as in Fig.^;. In the ID case, 
the thermal distribution leads to an even narrower line and 
higher rate, as small collision energies are enhanced by the 
density of states. 

In the quasi-2D or quasi-lD regime (Fig.^1), the line is 
almost unbroadened by temperature and recovers its natural 
Lorentzian shape in all theories. The temperature is so low 
that there is almost one collision energy in the system: 
in the 3D theory, ^hfl in the 2D theory, and HQ in the ID 
theory. 

In Fig. we consider a higher optical length (Z opt = 500 
Bohr), which can be reasonably attained using the intercom- 
bination lines close to thresold, or by increasing the photoas- 
sociation laser intensity. The line is now power-broadened, 
as predicted by Eqs. J46I47L while the free-space line is not. 
As explained before, the broadening is due to the zero-point 
momentum in the 2D case, and to fermionization in the ID 
case. In the quasi-2D regime, the line is also shifted by the 
logarithmic term, according to Eq. ( 145 \ . As suggested in 
Section Hi Gl this shift could be used as a way to measure 
the optical length. However, the line is progressively power- 
broadened as it is shifted. In Fig. [3] we plot the shift and 
width as a function of the optical length. Although the ex- 
periment may not be easy, this graph suggests that the shift 
should be observable. 

We now illustrate the theory with the intercombination 
lines in Strontium 86. The background scattering length of 
this species has been found to be quite large, between 610 
and 2300 Bohr |8J. As a result, the photoassociation line in 
the lattice is notably modified by the scattering length. In 
Fig. |4j we plotted the photoassociation rate for at, g = 600 
Bohr and a^g = 2000 Bohr. The effect of the scattering 
length is clearly noticeable, as predicted qualitatively by 
Eqs. J48I49I . As noted previousy, measurements of the rate 
for different lattice wave lengths might be able to give a bet- 




Detuning (kHz) Detuning (kHz) 

Figure 1: Photoassociation rate of Strontium 88 at low intensity (Z opt = 1 Bohr ~ 0.0529 nm) as a function of laser detuning, comparing 
free-space (3D), ID-confinement (2D geometry) and 2D-confinement (ID geometry), with a trapping frequency = 2-7T x 50 kHz. Black 
curves are the free-space 3D prediction - Eq. 13 1 i . Red curves are the effective 3D rate in a 2D geometry - Eq. 1251 . Green curves are 
the effective 3D rate in a ID geometry - Eq. 1291 . Panel (a) : T = 100 /iK (free-space regime), Panel (b) : T = 20 /j,K; Panel (c) : 
T = 4 fjK (confinement-dominated regime), Panel (d) : T = 0.1 /jK (Quasi-2D or - ID regime). In the confinement-dominated regime, 
the approximate line shapes calculated from Eqs. <34I36> (considering only the trap ground-state contribution) are represented by dashed 
curves. 



ter constraint on the scattering length. 

We also plotted on Fig. [4] the effective rate obtained if 
one replaces the effective scattering length j!5l > by the more 
usual resonant scattering length dl9l >. The latter appears to 
be insufficient to predict the correct rate for large scatter- 
ing length, showing that in this particular case, one cannot 
simply extrapolate the original work of Refs. Il4t[l5ll by re- 
placing the scattering length by the resonant scattering length 

IV. CONCLUSION 

We investigated how photoassociation is affected by the 
confinement of a ID or 2D optical lattice. We showed that 
effects are expected to be observable in the case of narrow 
resonances, like the intercombination lines found in alkaline- 
earth species. The main effects of confinement are: shift of 



the resonance (most importantly by the zero-point energy of 
the lattice), narrower lines, increase of power-broadening at 
low temperature, and direct sensitivity to a large value of the 
scattering length of the atoms. More specifically, a 2D con- 
finement could significantly increase the spectroscopic reso- 
lution, making it possible to observe the natural line width at 
relatively high temperatures. At low temperature, a ID con- 
finement induces a shift of the line which is proportional to 
the intensity of the photoassociation laser. This shift could 
be used to determine the strength of a resonance, without 
having to measure the density of the system. 
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Figure 2: Same as Fig. Qi for higher intensity or stronger resonance 
(/opt = 500 Bohr ~ 26.6 nm). The 2D and ID lines are now sat- 
urated (especially the ID line), in accordance with Eqs. 1461 and 
1471 . while the free-space line is not. The 2D line is also shifted 
according to Eq. 1451 . We indicated the three contributions to the 
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Figure 3: Shift AE of the photoassociation line in a 2D geom- 
etry with trap frequency Q = 2n x 50 kHz, as a function of 
optical length (proportional to photoassociation laser intensity). 
For convenience, we substracted the zero-point energy shift ^hQ, 
which is constant. The orange diamonds correspond to a tempera- 
ture T = 0.4 fiK, and the blue dots corresponds to a temperature 
T = 0.1 ^iK. These points are obtained by calculating the line shape 
using Eq. 1401 and finding the shift numerically. For each point, we 
represented a vertical line whose length corresponds to one sixth 
of the photoassociation line width; this gives an indication of how 
accurately the shift can be determined experimentally. The dashed 
lines are obtained using the approximate formula Eq. <45> . 
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